Intermittency and Local Heating in the Solar Wind 
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Evidence for inhomogeneous heating in the interplanetary plasma near current sheets dynamically 
generated by magnetohydrodynamic (MHD) turbulence is obtained using measurements from the 
ACE spacecraft. These coherent structures only constitute 19% of the data, but contribute 50% of 
the total plasma internal energy. Intermittent heating manifests as elevations in proton temperature 
near current sheets, resulting in regional heating and temperature enhancements extending over 
several hours. The number density of non-Gaussian structures is found to be proportional to the 
mean proton temperature and solar wind speed. These results suggest magnetofluid turbulence 
drives intermittent dissipation through a hierarchy of coherent structures, which collectively could 
be a significant source of coronal and solar wind heating. 
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The energy required to explain the coronal heating 
problem and, hence, the existence of a solar wind [l|, [2| 
might be provided by magnetohydrodynamic (MHD) tur- 
bulence [3, 4] driving a cascade [5] to scales where kinetic 
dissipation is efficient. This could account for the non- 
adiabatic expansion of heliospheric protons [6] . Analysis 
of higher order statistics from observations and numerical 
simulations indicates inertial range intermittency 0, is 
associated with coherent structures such as current sheets 
(oj-ITlj . These may play a role in heating the solar corona 
by magnetic reconnection (or nano flares). Connections 
between the turbulence cascade and intermittent dissi- 
pation are fundamental in hydrodynamics but are 
almost unexplored in space plasma physics, which relies 
heavily on linear Vlasov theory under the assumption of a 
uniform equilibrium plasma |13l4l5j . Recently, a statisti- 
cal link was found between coherent magnetic field struc- 
tures and elevated temperatures [l6| (also see llTl]), the 
interpretation of which has been questioned [isj. Here 
evidence is presented, based on correlations at two ranges 
of spatial scales and on global conditional averages, that 
temperature enhancements are locally linked to coherent 
structures and intermittency in solar wind turbulence. 

Motivation for this investigation lies at the core of tur- 
bulence theory. The Kolmogorov 1941 framework [19] 
provides a useful description, but a more precise theory 
must account for large fiuctuations in the energy dissipa- 
tion rate [20] . The resulting intermittency can be defined 
as the increasing non-Gaussian character of increments 
with decreasing scale [l2j. It is this connection between 
non-uniform dissipation and the statistics of increments 
that is embodied in the Kolmogorov Refined Similarity 
Hypothesis (KRSH) [21]. Since solar wind turbulence is 
well described by ideas that parallel its hydro dynamic 
antecedents — second order j22j , third order [23|, [2^ , 
and higher order statistics — it is perhaps paradoxi- 
cal that interplanetary dissipation is not presumed to be 
highly spatially non-uniform. Indeed, we are not aware 
of an alternative rationale that provides a physical expla- 
nation for solar wind intermittency. 



A major difficulty in linking energy dissipation to inter- 
mittency within low collisionality plasmas is that the un- 
derlying dissipation mechanisms are not unambiguously 
known. Often solar wind dissipation is framed in terms 
of Landau damping, cyclotron resonance, and instabili- 
ties in a uniform plasma where linear Vlasov theory is 
relevant. However, the observed consistency with inter- 
mittency theory implies deeper connections to KRSH. In 
order to pursue an unbiased assessment, without making 
assumptions about the dissipation function, proton tem- 
perature is adopted as a surrogate for the dissipation. 
Therefore, if the dissipation is non-uniform or intermit- 
tent, it will be refiected in the temperature distribution. 
In particular, we would expect to find sources of temper- 
ature enhancements within and around current sheets, 
which are the characteristic small scale coherent struc- 
tures in MHD turbulence [iliil,!!!. This is despite the 

collisional MHD framework not being applicable at the 
scales where the most intense current sheets would form. 
This lack of clarity regarding the balance between fiuid 
and kinetic phenomena is partly responsible for the wide 
disparity of viewpoints concerning coronal and solar wind 
dissipation. The use of proton temperature as a surrogate 
side-steps this issue and addresses the statistics of dissi- 
pation without detailed knowledge of the mechanism. 

We analyze the entire 64 s resolution magnetic field 
and proton temperature datasets from the MAG [26^ and 
SWEPAM [27] instruments onboard the ACE spacecraft 
at 1 AU. Rapid changes in the magnetic field are de- 
scribed by vector increments: 



Ab(t, At) = h(t + At) - b(t) 



(1) 



where b(t) is the magnetic field time series and At is the 
time lag. Here the temperature data cadence defines the 
lag At = 64 s employed. Using Taylor's hypothesis this 
corresponds to a spatial separation within the inertial 
range, which extends from about an hour to a second 
in the spacecraft frame [sj. In order to identify coherent 
(non-Gaussian) structures in the solar wind, a time series 
of the normalized partial variance of increments (PVI, 
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denoted I in Fraktur script) is constructed: 

lAbl 



(|Ab|^) 



(2) 



where (• • •) denotes an ensemble average. Events are 
selected by imposing thresholds on PVI, leading to a 
methodology that is comparable to classic magnetic dis- 
continuity identification in both solar wind observations 
and numerical simulations 1281. 
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FIG. 1. Mean Tp conditioned on the distance from (a) events 
identified using the PVI statistic. These events are linked 
to elevated temperatures. The temperature profile associated 
with < ^ < 1 is reproduced (b) to highlight its local mini- 
mum. These low value fiuct nations are heated upon encoun- 
tering current sheets. For reference, the dashed line represents 
the mean Tp for the entire dataset. 

In order to study heating near discontinuities, we com- 
pute averages of Tp conditioned on the distance from PVI 
events which exceed a threshold value 6: 

fp{Ar, 0) = {Tpir^ + Ar)\0 < ^{r^) < ^ 1) (3) 

where rc^ is the position of a PVI event and Ar is the 
spatial lag measured relative to rc^, obtained assuming 
frozen-in flow such that Ar ^ —Vs^At. 

Figure 1 shows the conditionally averaged proton tem- 
perature in the vicinity of a PVI event for selected values 
of 0. For thresholds above and including = 2, there is a 
peak in the average Tp at the location of discontinuities, 
as established in [16]. The steep fall-off at nearby spatial 
lags produces a deflned local maximum with a width of 
about 10^ km, which is equivalent to around a tenth of 



the turbulence correlation scale |2£ 



This local T„ 



maximum is more pronounced near stronger discontinu- 
ities, as seen when the PVI threshold is raised. The av- 
eraged Tp then transitions to a gently declining plateau, 
which is distinct for each value since the strongest 
PVI events are associated with the highest temperatures. 
Analysis of the data suggests that once a strong discon- 
tinuity is identifled, there is an elevated probability of 



flnding additional strong events within several correla- 
tion lengths (of order 10^ km). The net effect of this 
clustering, or non-Poisson property [see 30, for details], 
of strong PVI events is that the surrounding plasma is, 
on average, hotter. 

The cumulative mean waiting distances between non- 
Gaussian events are consistent with the interpretation of 
clustering PVI events. Beyond a spatial separation of 
around 0.1 Ac, there is an increased likelihood of encoun- 
tering another discontinuity. Therefore, we interpret the 
central peak of Tp enhancement in Fig. 1 as the result 
of local heating by individual coherent structures, while 
the broader plateau of elevated temperature is the result 
of an increased nearby density of strong coherent struc- 
tures. This clustering of heating events is indicative of 
a correlated inertial range intermittency process, in con- 
trast to an uncorrelated Poisson heating mechanism. 

It should be noted that all the conditional average Tp 
proflles with ^ > 1 lie well above the unconditioned aver- 
age value. Hence, it follows that the proflle correspond- 
ing to low value fluctuations (in the range < ^ < 1) 
must lie below the average Tp value, as shown in Fig. 1. 
However, it is not obvious that these fluctuations should 
produce a central local minimum in the conditional av- 
erage temperature proflle. These low level fluctuations 
are the smoothest regions and represent the closest ap- 
proximation to uniform plasma conditions [28| . Since the 
most uniform samples are cooler than the surrounding 
plasma, our results suggest the dominant sources of tur- 
bulence heating are unlikely to be found using methods 
that assume a uniform plasma. 

Figure 2 shows the conditionally averaged Tp depen- 
dence on the magnitude of spatial separation from the 
central discontinuity |Ar|, for three PVI thresholds. This 
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FIG. 2. Mean Tp conditioned on the distance from events 
identified using the PVI statistic. The plots consist of a core 
component of enhanced temperature and an extended heated 
outer component. These are both well approximated by linear 
scaling and the obtained temperature gradients are steepest 
near the strongest (highest 0) coherent structures. 
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expanded scale allows the two-tiered structure of the av- 
eraged temperatures to be more manifest. The data is 
well approximated by linear fits, and the obtained tem- 
perature gradients confirm the presence of an enhanced 
Tp inner core and an extended (> 0.1 Ac) heated region. 

In order to study the macroscopic behavior of inter- 
mittent heating, the data is divided into 10 hour samples 
and in each sample we compute the fraction of data rep- 
resented by strong events satisfying ^ > 4. Figure 3(a) 
shows PDFs of Tp conditioned on this density of coherent 
structures. As the percentage of strong PVI events in- 
creases, the probability density increases for the highest 
average Tp and vanishes for the lowest. This is consistent 
with coherent structures being directly associated with 
enhancements in temperature. The 10 hour samples are 
sorted by average solar wind speed into bins of width 
100 kms~^, and the average coherent structure density 
in each bin is plotted as a function of the correspond- 
ing average Tp. Figure 3(b) shows the higher average Tp 
samples also contain a greater density of strong coher- 
ent structures. A remarkable feature is that the data is 
ordered automatically into a sequence of increasing solar 
wind speed Vgw Faster wind has higher Tp and a higher 
density of discontinuities [18]. The bin corresponding to 
the lowest wind speeds {vqw < 300 kms~^) lies on the 
same trend line, but has a lower than average Tp. The 
implied mutual correlation between Tp, Vsw^ and heating 
rate is completely consistent with a turbulence-heated 
corona and solar wind [e.g. .31j. 
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FIG. 3. (a) PDFs of Tp conditioned on the density of coherent 
structures with ^ > 4. As the percentage of strong PVI events 
increases, the probability density increases for the highest av- 
erage Tp and vanishes for the lowest, (b) The higher average 
Tp samples contain a higher density of strong PVI events, and 
are ordered according to increasing average solar wind speed. 

The evidence presented here suggests solar wind 
plasma is hotter in the immediate vicinity of discontinu- 
ities. The sharpest discontinuities are linked to the steep- 
est temperature gradients and the clustering of disconti- 
nuities leads to regionally elevated temperatures. These 



conclusions contrast those of [l8| who employed a dif- 
ferent methodology, and only studied the nearest points 
neighboring the local maximum. Our results suggest re- 
lationships involving the dissipation averaged over a 
scale s, the increments (on the same scale) of the MHD 
Elsasser fields dzs^ and the local temperature Tp. Sym- 
bolically we consider two relationships: 

es--^^Tp (4) 

The first is an adaptation of the Kolmogorov Refined 
Similarity Hypothesis to MHD [e.g. 32], where ~ refers 
to statistical distributions that result in proportionality 
of moments after appropriate averaging. From KRSH, 
which is widely supported in hydrodynamics but un- 
proven [12], emerge properties such as multi-fractal scal- 
ing due to the irregular spatial distribution of e^. In 
astrophysical plasmas there is generally no direct mea- 
sure of eg, although the present study does demonstrate a 
statistical connection between increments and local tem- 
perature enhancements: the second relationship in EqlH 
Given the plausibility of a link between the local statis- 
tics of Tp and e^, our results support the hypothesis that 
MHD turbulence drives intermittent dissipation in the 
solar wind through kinetic processes operating within a 
hierarchy of intermittent coherent structures. 

It is instructive to determine if the observed inhomo- 
geneous elevations in Tp could signal an underlying heat 
source that contributes significantly to solar wind heat- 
ing. An initial assessment is provided by estimating the 
contribution from PVI events to the total internal energy 
U (X ripTp^ where Up is proton number density. Here the 
PVI events are conditioned on threshold ^ > ^. The in- 
ternal energy U{0) and volume V{0) are associated with 
the central PVI event and the two nearest and next- 
nearest neighbors, which together form a local maximum 
(see Fig. 1). Figure 4 plots the conditional internal en- 
ergy density Ud{0) = U {6)/V{0) normalized to the inter- 
nal energy density of the entire dataset Ud{0 = 0). As 
the PVI threshold is raised, Ud{0) increases as is con- 
sistent with internal energy being concentrated in strong 
discontinuities. In order to better quantify the contribu- 
tion of coherent structures to we define non-Gaussian 
events as those where the ratio of PVI to Gaussian distri- 
butions exceeds an order of magnitude. For these events, 
the PVI threshold dependence of U and V is examined. 
Figure 4 shows events where ^ > 2.4 make up 50% of the 
total internal energy but only occupy 19% of the volume. 
The strongest PVI events, ^ > 5, occupy just 2% of the 
volume but contribute 11% of the internal energy. In a 
separate study [33], the threshold = 6 was found to 
identify events that are likely strong active magnetic re- 
connection sites. Therefore, it seems the inhomogeneous 
features identified here could contribute substantially to 
the internal energy budget of the solar wind. 

If the dissipation function is intermittent, as our re- 
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FIG. 4. (a) PDFs of PVI statistic from solar wind dataset 
compared to a unit- variance Gaussian distribution, where 
above = 2A (dot-dash hne) corresponds to non-Gaussian 
coherent structures, (b) Estimates of the fraction of total in- 
ternal energy density Ud = U/V in structures with ^ > ^ 
are found to increase with 0. This suggests internal energy is 
concentrated in discontinuities. 



suits suggest, its properties can be better understood in 
the context of large scale turbulence, which may pro- 
vide the framework to understand coronal and solar wind 
heating 0,0, [HI. Hence, variations in intermittent dis- 
sipation might account for differences in wind speed [31]. 
If demonstrated, this would bring closure to the inter- 
relationships between speed, temperature, and density 
of coherent structures, that remain otherwise presently 
unexplained. As observations probe below the ion in- 
ertial scale to study the kinetic physics of dissipation 
[14 . 36] and closer to the solar wind sources in the Solar 
Probe and Solar Orbiter missions, the link between iner- 
tial range increments and non-uniformity in dissipative 
structures will become more important. Indeed, such ef- 
forts might establish a much closer relationship between 
solar wind turbulence and the classical hydrodynamic 
concept of intermittent dissipation [20, 21]. 
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